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ABSTRACT 



^ \ The HST archival UV imaging polarimetry data of NGC 1068 is re-examined. 

t;;-!- \ Through an extensive estimation of the observational errors, we discuss whether the 

cn ' distribution of the position angles (PAs) of polarization is simply centrosymmetric or 

^ . not. Taking into account the effect of a bad focus at the time of the observation, we 

Q\ \ conclude that, within the accuracy of HST/FOC polarimetry, the PA distribution is 

^ I completely centrosymmetric. This means that the UV polarization originates only from 

,-^ ■ scattering of the radiation from a central point-like source. 

However, our analysis shows that the most probable location of the nucleus is only 

O I ~ 0."08 (~ 6 pc) south from the brightest cloud called "cloud B". The error circle of 

^ I 99% confidence level extends to cloud B and to "cloud A" which is about 0."2 south of 

cd ■ cloud B. By this FOC observation. Cloud B is only marginally rejected as the nucleus. 
Assuming that the UV flux is dominated by electron-scattered light, we have also 

k/< I derived a three-dimensional structure of the nuclear region. The inferred distribution 

;h ' suggests a linear structure which could be related to the radio jet. 

Subject headings: galaxies:Seyfert, galaxies:individual (NGC1068), ISMxlouds, meth- 
ods:data analysis, polarization, scattering 



^ Based on observations made with the NASA/ESA Hubble Space Telescope, obtained from the data Archive at the Space Telescope 
Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 
5-26555. 



1. Introduction 

The nucleus of the Seyfert 2 galaxy NGC 1068 
is now firmly believed to be obscured from direct 
view. This was clearly inferred by Antonucci & Miller 
(1985) from a spectropolarimetric study, and now one 
of the major problems in this galaxy is to locate this 
hidden nucleus very accurately on the high-resolution 
images taken by the HST and on VLBI radio maps. 
This is particularly important, since the investigation 
of the physical conditions and kinematics of the nu- 
clear vicinity is greatly influenced by the exact loca- 
tion of the hidden nucleus. 

Several authors have addressed this problem, but 
most of the nuclear positions determined are indirect 
ones in the sense that they are from lower resolu- 
tion images than the ffST images or VLBI maps, and 
these positions are slightly different from one another 
(see e.g. Thatte et al. 1997 and references therein). 
The only exception so far for pinpointing the loca- 
tion directly on the HST high-resolution images is to 
use imaging polarimetry data. The nucleus can be 
determined as the center of the centrosymmetric dis- 
tribution of the position angle (PA) of polarization, 
which is expected to be observed if the radiation from 
the nucleus is being scattered by the surrounding gas. 
(We refer to this case as a 'point-source scattering 
case' hereafter.) Capetti et al. (1995a, b) have de- 
termined the nuclear location by this method. How- 
ever, in their PA map, clear deviations from the cen- 
trosymmetric pattern are seen, and they did not state 
whether they are real, or discuss the observational er- 
ror in PA at each position of the image. If the de- 
viations are real, the nuclear position determined by 
them would not be very accurate and reliable. 

In this paper, we extend their work to discuss these 
deviations through the extensive estimation of the ob- 
servational errors, and then to re-determine the loca- 
tion of the nucleus more accurately within a convinc- 
ing error circle. Our result is that the most probable 
location of the nucleus has moved to the north by 
~ 0."2, significantly larger than the quoted error of 
Capetti et al. (1995b). 

The second objective of this paper is to derive 
the three-dimensional structure of the nuclear region. 
The HST images show that this region has a knotty 
and filamentary structure. If the UV radiation from 
each knot is dominated by the scattered radiation, the 
polarization degrees of the clouds provide the scat- 
tering angle at each cloud. Hence, we can derive the 




three-dimensional distribution of these clouds with re- 
spect to the nucleus. 

We describe the data in § || and error estimation 
In § 0, we examine the PA distribution, and 
5[ we discuss the location of the nucleus. Then 
6l we derive the three-dimensional distribution 
of the scatterers. We discuss these results in § M and 
our conclusions are presented in § R In the appendix, 
we summarize the method of error estimation. We 
assume a distance of 14.4 Mpc to NGC 1068 in this 
paper (TuUy 1988), corresponding to a scale of 1" ~ 
70 pc. 

2. The archival data 

2.1. Reduction procedure 

The archival HST data used are summarized in Ta- 
ble 0. The polarimetry data were obtained on Febru- 
ary 28, 1995. These data have been pubhshed by 
Capetti et al. (1995b). The central small portion was 
in the 10%-level nonlinear regime, so a flat-field lin- 
earity correction was applied, but we confirmed that 
this correction has no significant effect on our anal- 
ysis, by implementing the same analysis below on 
the data without a linearity correction. The data 
were subsequently processed in the standard manner 
to correct for geometric distortion and flat-field re- 
sponse. The reseau marks were removed using neigh- 
boring pixels. The background subtraction was im- 
plemented using the outermost regions of the images. 

These data were obtained after the COSTAR de- 
ployment, but at the time of this polarimetry observa- 
tion HST had an extraordinarily poor focus due to a 
large movement of the secondary mirror, the effect of 
which will be discussed in detail below. Therefore, the 
images through different polarizers (POLO, POL60, 
POL120), which are known to be slightly shifted rel- 
ative to one another, were registered by using the im- 
age shift values of calibrations by Hodge (1993, 1995), 
not by using point-like sources in the outer region. 
Then the images were scaled according to the expo- 
sure times and the transmittances of the polarizers 
and the F253M filter. (The filter transmits the UV 
radiation around 2400-2700 A.) Finally these images 
were combined to obtain the Stokes parameter /, Q, 
U images. 

Using several archival HST/FOS spectra taken 
with 0."3 aperture, we have estimated the amount 
of the emission-line contamination in the F253M im- 
age to be around the 10 ^ 15% level. This includes. 



Table 1 
hst/foc archival data used 



Rootnamo 


Filter 1 


Filter 2 


Exp. time (sec) 


Date of obs. 


Description 


x274020at 


F253M 


POLO 


1796.625 


Feb 28, 1995 


UV continuum 


x274020bt 


F253M 


POLO 


341.625 


Feb 28, 1995 




x274020ct 


F253M 


POLO 


1451.625 


Feb 28, 1995 




x274020dt 


F253M 


POL60 


748.625 


Feb 28, 1995 




x274020et 


F253M 


POL60 


1044.625 


Feb 28, 1995 




x274020ft 


F253M 


POL60 


1201.625 


Feb 28, 1995 




x274020gt 


F253M 


POL60 


591.625 


Feb 28, 1995 




x274020ht 


F253M 


POL120 


1608.625 


Feb 28, 1995 




x274020it 


F253M 


POL120 


1796.625 


Feb 28, 1995 




x24e0102t 


F501N 




1196.000 


Jan 10, 1994 


[OIII] 


x24e0103t 


F501N 


F4ND 


1196.000 


Jan 10, 1994 




x24c0106t 


F253M 




1196.000 


Jan 10, 1994 


UV continuum 



however, the broad Fell Hnes which are also a part of 
the scattered hght in addition to the scattered con- 
tinuum (Antonucci, Hurt, & Miller 1994). Therefore, 
the effect of the line contamination on the polariza- 
tion should be much smaller. 

We have also used the archival images of NGC 1068 
taken through the F501N filter, i.e., the [OIII] image, 
and through the F253M filter but without polarizers, 
the focus of which seem to be fine. They were ob- 
tained on January 10, 1994, also after the COSTAR 
deployment. They have been published by Macchetto 
et al.(1994). The images were processed through the 
same procedure as above, but in the F501N image, 
the central nonlinear and saturated portion was filled 
with the appropriately scaled image taken with the 
F4ND filter. The pattern noise associated with the 
FOG nonlinearity was removed by Fourier filtering. 
The registration of these two images was carried out 
using point-like sources in the outer region, with the 
uncertainty of less than 1.4 pixel. The registration of 
this F253M image and the / (total intensity) image 
from the polarimetry data was done by taking a cross- 
correlation, because of the bad focus in the / image. 
We also tried this registration using point sources and 
found the results coincide within 1 pixel. 

Both of the observations were implemented in nor- 
mal 512 X 512 mode, where the pixel size is 0."014 x 
0."014 and the field of view is 7" x 7". 



2.2. HST Focus and Degradation of the Image 

Due to the bad focus in the imaging polarimetry 
observation, the images have a significant blurring. 
The effect on the polarization analysis should be sig- 
nificant especially at the region where the gradient of 
the polarized flux distribution is large. We first esti- 
mate the extent of this blurring by using the fine- focus 
(non-polarimetric) F253M image. 

To avoid confusion, we call this F253M image with- 
out polarizers a "sharp" UV image, and the / image 
of the polarimetry simply an / image. Synthetic aper- 
ture photometry was carried out to take the ratio of 
the / image to the sharp UV image at various aper- 
ture diameters, and the results on a few tens of local 
maxima over the image were averaged. The results 
are shown in Figure m The scatter of the ratio for 
each aperture is shown as a vertical line. The ex- 
pected amount of the error in the ratio due to statis- 
tical noise is also shown with horizontal ticks. As we 
expect that the local peak intensity is lowered by the 
blurring, the average ratio is certainly less than unity 
and it becomes smaller with smaller aperture, about 
0.9 at 0."15 aperture, down to ~0.85 at 0."05 (though 
the statistical significance of these average ratios be- 
comes lower with smaller aperture). The regions in 
the / image with this ratio larger than ~ 1.0 have sig- 
nificant 'leaks' from the neighboring bright regions. 

To infer correctly the degradation of the point 
spread function (PSF) of this observation, we have 
estimated the amount of the secondary mirror motion 
in the following way. We constructed model PSFs us- 
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Fig. 1. — The averaged ratio of the I image to the 
sharp UV image over various aperture diameter is shown 
as squares, connected by dotted lines. The vertical line 
on each square represents the standard deviation of the 
scatter from each averaged ratio. The expected amount 
of the error in the ratio due to statistical noise is shown 
in horizontal ticks. 

ing TinyTim (Krist 1993) with a range of secondary 
mirror movement and convolved the sharp UV image 
with those PSFs. Then we compared them to the / 
image which is convolved with the PSF of nominal fo- 
cus for equivalence. We calculated the ratio of these 
two by the same synthetic aperture photometry as in 
Figure |[ It is expected that this ratio becomes unity 
for all aperture sizes when the appropriate amount of 
the secondary mirror shift is taken for the model PSF. 
We found that this was the case for a shift of about 
9 ~ 10 microns. The model PSF at this focus indi- 
cates that the FWHM of the PSF was 0."10 - 0."15 
at the time of the observation. 

Therefore, we decided to implement our analysis in 
three ways; (A) 10 x 10 pixel binning, corresponding 
to ~ 0."15 X 0."15, (B) 20 X 20 pixel binning, (C) 
0."15 aperture diameter synthetic photometry. The 
size for the cases (A) and (C) could be too small, but 
the effect of small aperture will be taken into account 
in the error estimation. For case (C), the aperture 
centers are selected as the local peaks in the / image, 
which are thought to represent the positions of the 
resolved clouds in this nuclear region. 



3. Error estimation 

In order to examine the PA distribution accurately, 
we need extensive estimation of the observational er- 
ror in the HST/FOG polarimetry. We describe major 
ideas here and summarize other details in appendix. 

We estimate the error in polarization degree and 
PA by considering the following four sources of errors 
that are expected to be major among various error 
sources. 

(1) Statistical error. Poisson noise is assumed. 

(2) Uncertainty from the image reg istr ations of 
three polarizer images. As described in 



2.1 



we have 

registered the images using the calibrated values of 
Hodge (1995) and their wavelength dependence by 
Hodge (1993). The calibrated image shifts have an 
uncertainty of about ±0.3 pixel. Therefore we shifted 
POL60 and POL120 images by this amount along the 
image x and y axes relative to POLO image and cal- 
culated the resulting change of the polarization. 

(3) Polarizer axes uncertainty. The transmission 
axes of each polarizer have nominal values of 180°, 
60°, 120° with respect to the +x axis of the image, 
and those values have about ±3° uncertainties (Nota 
et al. 1996). The errors in the calculated polarization 
due to this uncertainty were estimated by considering 
the dependence of the calculated polarization on the 
polarizer axes' angles (see appendix for more detail). 

(4) We consider the uncertainties from the follow- 
ing error sources as the uncertainties in the "correc- 
tion factors" to be multiplied to each image through 
each polarizer. 

(i) Each polarizer, especially POL60, has differ- 
ent wavelength dependence of transmittance, so 
when the spectrum of the object changes from 
place to place, each polarizer's image has to be 
multiplied by a different factor to correct for 
the different effective throughput. By assuming 
a power-law source spectrum and varying the 
power-law index in a reasonable range, we esti- 
mate that this has effects of less than about 1% 
in the correction factors. 

(ii) Differences in PSFs through each polarizer re- 
sult in uncertainties in these multiplying factors 
locally. We estimate this effect from the result 
given by Hodge (1995), who has investigated 
the differences in the fluxes of unpolarized point 
sources through different polarizers over various 



aperture sizes. Although this is for data with 
nominal focus, the maximum uncertainties for 
our out-of-focus data can be estimated for each 
of our binning size (A), (B) and (C) above, from 
the flux differences of smaller aperture size in 
Hodge's result. For cases (A) and (C), the rel- 
ative discrepancies between the fluxes through 
each polarizer are expected to be less than 3 to 
4 %, reading from Hodge's plots for 3 to 5 pixel 
aperture radius, and for case (B) less than 2 to 
3%, from the same plots for ~ 7 to 10 pixel 
aperture radius. The fact that Hodge's result 
is for point sources while our sources are more 
extended, gives us additional support that the 
values above can be considered as the maximum 
uncertainties. 

(iii) Flat-fielding uncertainty. The flat fields used 
by the FOC standard calibration are heavily 
smoothed and they do not correct for the fine 
scale features in the original, unsmoothed flat 
fields. We have estimated the resulting flat- 
fielding uncertainties by calculating the devi- 
ations of the unsmoothed fiat field from the 
smoothed one, for the used portion of the im- 
age. The results are about 3% for the cases (A) 
and (C), 2% for the case (B). 

Taking the sum of the squares of these (i)'^(iii) 
factors, we consider uncertainties in each multiply- 
ing factor as 5% for cases (A) and (C), 4% for case 
(B), at most. The effect of these uncertainties on the 
polarization were calculated from the dependence of 
the polarization on these multiplying factors (see ap- 
pendix for more detail). 

The total error was estimated by combining in 
quadrature all the uncertainties above from (1) to 
(4). Source (4) is usually larger than (2) and (3). 
Typically, the total estimated error in the polariza- 
tion degree, crp, in the regions of sufficient statistical 
S/N (the portion of the data used in the analysis be- 
low) was about 4.5 ~ 5.5% for the cases (A) and (C), 
while about 3.5 ^ 4.5% for the case (B). If the statis- 
tical error is excluded, the quadratic sum of the other 
three was about 4.5% for cases (A) and (C), about 
3.5% for case (B). 



4. Examination of Position Angle Distribu- 
tion 

As shown in Capetti et al. (1995b), the overall 
PA distribution is quite centrosymmetric, quite close 
to a 'point-source scattering case'. However, clear 
deviations from the centrosymmetric pattern are seen 
in some regions, especially around the very center (see 
their Fig. 2). We now examine these deviations and 
discuss whether the data are really consistent with 
this simple model, based on the PA errors estimated 
in the previous section and also on the investigation 
of the bad-focus effect discussed in § 2.2. 



4.1. Construction of x^ image 

We evaluate whether the data are really consis- 
tent with point-source scattering by calculating the 
X^ value. For each iih binned pixel or aperture, de- 
note the PA data as 6'pa(«), and the ideal centrosym- 
metric PA with a certain symmetric center (xc, yc) as 
6'idoai(i ;xc,yc)- We write the x^ value for (xcVc) as 

valid data r„ , ... „ ,. .-.o 

2, ^ _ Y^ {t'pA(«) - t'idcaU* ;a;c,yc)r 



X^(2^c,yc)= X! 



o-epA(«)^ 



(1) 

using the total estimated error in PA, crgp^^, described 
in the previous section. The x^ values for various 
symmetric centers are calculated and a 'x^ image' is 
constructed, in which the image value at a certain 
point represents the x^ value with this point being 
the symmetric center. 

The ideal centrosymmetric PA for one large binned 
pixel or aperture depends on, in addition to the direc- 
tion to the symmetric center, the polarized fiux distri- 
bution within that pixel. Therefore we used the orig- 
inal, not binned, distribution of polarized fiux within 
the binned pixel or aperture for the calculation of 
the x^ value. That is, we calculate the "polarized- 
flux-weighted" ideal centrosymmetric PA (actually, 
weighted further by the reciprocal of statistical error 
of polarized fiux, to avoid having too much weight on 
the pixels with low S/N polarized fiux) as 

6'ideai(i ; a;c, Vc) = 7:, arctan ' ^"^ , .. , (2) 



2 Qidcal(i) 



where 



;idcal 



(0 



E 



ith bin 



klp{j,k) cos2ejk 



Eith bin 



, (3) 




of pixel {j,k) 



\, symmetric center 



Fig. 2. — Illustrates how the centrosymmetric PA for 
each binned pixel is calculated using the polarized flux 
distribution within the binned pbcel as weight. Thick solid 
vectors represent the PAs perpendicular to the direction 
to the symmetric center. 



UidcaS) 
Wjk 
Ojk 



Effc '" Wjkip (j, k) sin 29 jk 



Eiih bin 



af^'{],ky 



arctan 



Vk -Vc 



, (4) 

(5) 
(6) 



Tiie quantities Ip{j, k), crf^^ij, fc), {xj,yk) are the po- 
larized flux, statistical error in the polarized flux, and 
central point, for the original pixel (j, k) within the 
ith binned pixel, respectively. These are illustrated 
in Figure g. The summations are taken over all pix- 
els {j, k) within the ith binned pixel. The angle Ojk 
corresponds to the ideal centrosymmetric PA for the 
center of the pixel (j, k). 

This correction in the ideal PA (the difference be- 
tween the polarized-flux-weighted ideal PA and the 
centrosymmetric PA which is calculated simply for 
the center of the binned pixel or aperture) was found 
to be very important in the central bright and knotty 
region. In this region, the amount of the correction 



was larger than the total error of PA estimated in § |[ 

4.2. The effect of bad focus on tfie evaluation 

ofx' 

The extraordinarily degraded PSF of these po- 
larimetry data can affect the PA distribution signifi- 
cantly, especially in the pixels surrounding the bright 
clouds, even with relatively large pixel binning. We 
eliminate this effect from the analysis by masking out 
the binned pixels or apertures which are suspected 
to have significant 'leak' from the surrounding pixels. 
If we take the ratio of the / image to the sharp UV 
image for each bin or aperture, the pixels or aper- 
tures which are not affected by the surrounding re- 
gions should have this ratio less than around unity, 
which can also be down to ~ 0.9 for 0."15 aperture 
as shown in Figure |^ and discussed in § 2^ . We call 
this ratio here a "contamination measure" . We should 
mask out the region where this contamination mea- 
sure is much larger than unity, so we set an appropri- 
ate threshold on this contamination measure, which 
we denote as [CMJ^^^^^. In addition, we use only the 
region where the polarization has been detected with 
high statistical S/N. We set an appropriate threshold 
on P/ap'^*' (statistical S/N of P), which we write as 

The result is that, choosing appropriate thresholds 
on these two factors, we have obtained minimum re- 
duced x^ (the X^ divided by degrees of freedom) of 
slightly less than unity for all three binning cases from 
(A) to (C). This indeed means that these polarime- 
try data are totally consistent with point-source scat- 
tering, within the accuracy of the FOC polarimetry. 
In Figure 0, we show the map of the minimum re- 
duced x^ value in the case (A), for various sets of 
PM]max and [S/N]^;^. The x^ becomes smaller for 
larger [S/NJ^^^j^^ and smaller [CMJ^^^^, and we clearly 
see the x^ becomes minimum and almost constant for 
[S/N]^i„ > 8 and [CM]„^^^ < 1.1, with its value indi- 
cating the fit is very good. The dotted-line contours 
indicate the cut-off value of the reduced x^ with 99% 
confidence level for each degree of freedom in each set 
of [S/N]j^ij^ and [CMJ^^^^^. The x^ also becomes stable 
for [S/N]^i„ > 8 and [CM]^^^ > 1.3, but the reduced 
X^ value for this region is rather hard to accept, since 
it is much larger than the cut-off value. Furthermore, 
there is a clear transition between these two regions 
where the gradient of x^ is large and x^ switches to 
an acceptable value when [CMJ^^^^^ becomes smaller. 
Similar results have been obtained for the cases (B) 
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Fig. 3. — The map of the minimum reduced y^ for case 
(A) is drawn in solid contours on the plane of the thresh- 
olds on the contamination measure (ratio of the I image 
to the sharp UV image) and statistical S/N (Pjap). Also 
shown in dotted contours are the cut-off values of the re- 
duced x^ with 99% confidence level for each degree of 
freedom in each set of thresholds. 

and (C). Therefore we conclude that the PA distribu- 
tion is consistent with a point-source scattering only 
when the contaminated regions are excluded. 

5. The location of the nucleus 

Now that the point-source scattering fit is proved 
to be acceptable, we can discuss the most probable 
location of the nucleus as the position of the minimum 
X^ point, and its error as the contour of the x^ image. 

5.1. Minimum x^ point 

We show the PA distribution and the point with 
the x^ being minimum for each case (A), (B), and (C) 
in Figure 0, H, and 0, respectively. The observed PAs 
for each binned pixel or aperture are shown in two 
white vectors that correspond to 6'pa ± crgp^ , where 
agp^ is the total estimated error in PA. Also shown 
in black vectors are the polarized-flux-weighted ideal 
centrosymmetric PA with the symmetric center being 
the minimum x^ point, marked as a plus, for each bin- 
ning case. For these figures, we have chosen [CMJ^^^^^ 
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Fig. 4. — The PA distribution with 10 x 10 pixel binning 
[case (A)] is shown on the / image of the UV polarimetry. 
The position of minimum x is shown as a plus sign, and 
the error circle of 99% confidence level is also drawn. The 
two white vectors for each binned pixel indicate ^pAicrgp^ 
and the black vector shows polarized-fiux-weighted ideal 
centrosymmetric PA with the minimum x^ position being 
the symmetric center. The image is in log scale, and has 
been cut at 0.01% of the peak intensity. 

to be 1 in all cases, which is considered to be a natu- 
ral constraint on the contamination measure (though 
we can also take the values of up to '--^ 1.1 to allow 
for the statistical noise, as seen in Figure 0). We have 
taken the lowest [S/Nj^^j^^ in each binning case that 
could yield acceptable y^ values, in order to include as 
large an area as possible. The minimum reduced x^ 
value was found to be 0.88 for 41 degrees of freedom 
with [S/N]^i„ = 8 in case (A), 0.80 for 26 degrees with 
[S/N]j„in = 7 in case (B), and 0.74 for 28 degrees in 
case (C). The underlying image in Figure Q and Fig- 
ure g is the / image, while the sharp UV image is 
presented in Figure g for comparison. 

Compared to the previous result from the same 
data (Capetti et al. 1995b), the minimum x^ point, 
i.e., the most probable location of the nucleus, is very 
close to the UV brightest cloud, in all cases from (A) 
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Fig. 5. — The same as FigW but with 20 pixel binning 
[case (B)], and the underlying image is the "sharp UV 
image", which was taken when the HST focus was fine. 
The image is in log scale, and has been cut at 0.01% of 
the peak intensity. 

to (C), as seen in Figure covering just the central 
'^ 0."5 region. The white dash-dot contour represents 
the UV image, while underlying is the [OIII] image 
through the F501N filter. Our minimum y^ points 
shown in small plus signs are only about 0."12 south 
from the UV peak at the origin of the coordinates, 
whereas Capetti et al. (1995b) have located the nu- 
cleus at 0."3 south from the UV peak (the value is 
actually taken from Capetti, Macchetto, & Lattanzi 
1997), shown as a large plus sign with its size being 
their estimated error. 

The reason for this difference is quite clear. Corre- 
sponding to the "transition" of the minimum x^ value 
in Figure Ijfrom [CMJ^^^^ of 1.3 to 1.1, we have found 
that the position of the minimum x^ point experi- 
ences the transition from ^ (0."01, — 0."23), which is 
almost within the error of the Capetti et al. point, 
to our points - (-0."03, -0."12). We found two re- 
gions on the plane of ([CMJ^^^^^, [S/Nj^^j^^) where the 
position of the minimum y^ point is stable, just cor- 
responding to the two regions in Figure ^ where the 



Fig. 6. — The same as Fig.kl but for the case of aper- 
ture photometry on the local peaks with 0."l5 diameter 
[case (C)]. The underlying image is the I image from the 
imaging polarimetry. It is in log scale, and has been cut 
at 0.1% of the peak intensity to show the bright regions 
more clearly. 



minimum x^ values are almost constant. These re- 
sults are for case (A), but we see quite similar results 
for cases (B) and (C). Therefore, the difference be- 
tween the two locations is expected to arise from the 
masking-out procedure for the regions contaminated 
by the focus effect. 

The reason for the shift of this direction is also 
fairly clear. We show the PA distribution in the 
masked-out region in Figure R enclosed by white 
lines, in addition to the PA distribution already shown 
in Figure H. The masked-out region, especially around 
the UV peak, seems to direct the nucleus to the south. 
So the minimum x^ point moves to the north when 
these regions are excluded. 

These PAs in the masked-out regions clearly can- 
not be fit by point-source scattering, and we have 
certainly shown that the minimum y^ becomes ac- 
ceptable only if we mask out these regions. The ten- 
dency is that the deviations in these regions can be 
explained by the contamination from the neighbor- 
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Fig. 7. — The enlarged view of the central ^ 0."5 region. 
The positions for minimum x^ points and 99% confidence 
level contours are shown in white dashed, solid, and dot- 
ted line for cases (A), (B), and (C), respectively. The 
black solid contour represents 99% confidence level error 
circle for case (B) with smaller uncertainty in the "cor- 
rection factor". Underlying is the F501N [OIII] image, 
and the dash-dot contours represent the sharp UV image. 
Both are in log scale. The thin solid and dotted plus sign 
indicates the nucleus location suggested by Capetti et al. 
(1995b) and Braatz et al. (1993), respectively, with their 
sizes being the errors estimated by the respective authors. 

ing bright regions. Since our niasking-out procedure 
is quite natural and reasonable, and this tendency 
gives us further support for our procedure, we con- 
clude that the Capetti et al. point should be revised 
to a more northern point, indicated by our minimum 
X^ positions. 

5.2. Error circle 

The error of the location of the symmetric center 
can be estimated as contours of the x^ image. The 
contour of 99% confidence level has been drawn for 
each binning method in Figures U ^ |6[ and these three 
different contours are all shown with an enlarged scale 
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Fig. 8. — The same as Fig. H but the observed PAs in 
the masked out region, enclosed by white lines, are also 
drawn. The small plus sign is the same as in Fig. M, while 
the large plus sign represents the location of the nucleus 
suggested by Capetti et al. (1995b), with its size being 
their estimated error. 

in Figure ^ in three different white circles. These 
three are slightly different from one another. This is 
partly due to the uncertainties in our estimation of PA 
errors, but the error circle for case (B) is smaller in 
the north-south direction partly due to the fact that 
the valid data points in case (B) are more extended 
to the east and west (Fig. a), while in cases (A) and 
(C) they are only extended to the north and south. 

As we described in § 0, the uncertainty in the "cor- 
rection factor" should be considered as the maximum 
possible, which might be indicated also by the fact 
that our minimum x^ value is slightly smaller than 
unity. We have also implemented the calculation of 
X^ with a slightly smaller value for this uncertainty, 
and found that the minimum x^ points and northern 
part of the error circles were almost stable, whereas 
the southern part of the error circles moved slightly 
to the north. This shift is about 0."05 if we take the 
uncertainty in the correction factor to be 4% for cases 
(A) and (C), and 3% for case (B), instead of 5% for 



Table 2 
Positions of the Nucleus and Other Sources 



source name 



offset west ("Y offset north {"Y 



R.A. (J2000) 



Dec. {J2000) 



UV and optical peak 
hidden nucleus 
SI source 



0.00 

-0.03 ±0.05 
0.02 ±0.10 



0.00 

-0 12+°°'^ 
-0.13 ±0.10*^ 



02'i42™40.=711±0.=005 
02'^42™40.=713±0.=006 
02h42™40.=710±0.=001 



-00°00'47."81 ± 0."08 



-00°00'47."93 



-o."ii 



^+0."14 

-00°00'47."94±0."02 <* 



''Offset from UV and optical peak. 

''From Capetti et al. 1997. 

'^Quadratic sums were taken to obtain errors. 

''From Muxlow ot al. 1996. 



(A) and (C), and 4% for (B). In Figure the error 
circle of 99% confidence level has been drawn as a 
black contour for the case (B) with 3% uncertainty in 
the correction factor. The minimum x^ values were 
1.16, 1.06, and 0.99 for the cases (A), (B), and (C), 
respectively, which are still low enough to accept the 
fit. Therefore, the error circle would be smaller in the 
southern part. 

The brightest [OIII] cloud in Figure is called 
"cloud B" and the fainter cloud 0."2 to the south is 
called "cloud A" (Evans et al. 1991; see also Bland- 
Hawthorne et al. 1997). Cloud B is slightly but signif- 
icantly displaced from the brightest UV cloud, where 

the offset is about 0."05 at PA 37°. These are 

consistent with the values noted by Macchetto et al. 
(1994). Note that the accuracy of the registration 
of these two images is estimated to be better than 
1.4 pixel, corresponding to 0."020, as described in 
The error circles suggest that the nucleus is 
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located between clouds A and B. The projected dis- 
tance between cloud B and our minimum x^ point is 
only ~ 0."08 (^ 6 pc). Cloud B is only marginally 
rejected as the symmetric center. This means both 
that the quality of these FOC polarimetry data can 
only limit the location with this amount of error, and 
that the nucleus could be located just beside cloud B. 
For comparison, the location of the nucleus de- 
termined by Braatz et al. (1993) as the peak po- 
sition at 12.4/im has been drawn in Figure 0, with 
the size being their estimated error. This location 
is given with respect to the optical continuum peak, 
but we have used their value with respect to the UV 
peak instead, since the UV peak and optical peak 



are coincident within the accuracy of the registra- 
tion of HST/WFPC2 and FOC images (Capetti et 
al. 1997; Kishimoto 1998). Our minimum x^ posi- 
tions are marginally within their error box. Thatte et 
al. (1997) have also located the nucleus as the peak 
of near-infrared emission, which is slightly north of, 
but almost identical to the position of Braatz et al. 

In conclusion, allowing for the displacement of the 
minimum x^ points and error circles of the four cases 
in Figure 0, we locate the nucleus 0."12 south and 
0."03 east from the UV peak, with the error circle 
extending from this point about 0."07 to the north, 
0."12 to the south, and 0."05 to the east and west. 
In Table 0, we summarize this result with equatorial 
coordinates, using the result of absolute astrometry 
by Capetti et al. (1997) on the fiST continuum peak. 

6. Three-dimensional structure 

We have shown that the UV polarization is re- 
ally consistent with point-source scattering and deter- 
mined the location of the nucleus accurately. Based 
on these two results, we infer the three-dimensional 
structure of the nuclear region, by using the polariza- 
tion degrees as the indicator of scattering angles or 
'viewing angles' at each resolved cloud in the image. 

6.1. Assumptions 

Interpretation of the polarization degree is usually 
difhcult, due to the existence of unpolarized, "di- 
luting" radiation in terms of polarization, such as 
starlight. In the nuclear region of NGC 1068, how- 
ever, there is some evidence that suggests the diluting 
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Fig. 9. — Illustrates the configuration of the cloud model, 
(a) Definition of viewing angle of a spherical cloud, (b) 
The two 'critical' positions determined from the polariza- 
tion of the cloud. 

radiation in the UV range is fairly small, as described 
below. If this is correct, we are in the unusual situa- 
tion that we can infer the distribution of the nuclear 
resolved clouds three-dimensionally. The polarization 
degree reflects the angle between our line of sight to 
the scatterers and the line of sight of those scatter- 
ers themselves to the illuminating source, which we 
call here 'viewing angles' of the scattering clouds, as 
shown in Figure ^. If we obtain the viewing angles 
of each cloud, we can locate each cloud along the line 
of sight with respect to the central radiation source, 
and thus we will be able to have a three-dimensional 
view of this nuclear region. 

We simply assume here that the UV radiation 
consists of only radiation scattered by free electrons 
and no other diluting component exists. This is 



based on the following two observational arguments 
in NGC 1068. 

(i) HST/FOS UV spectropolarimetry shows that 
the polarization is almost constant over the UV range 
and its degree declines only in the wavelength longer 
than -- 2800A (Antonucci et al. 1994). This has been 
interpreted to mean that the starlight contribution 
is significant only in the redward of 2800A, almost 
outside the F253M filter transmission. 

(ii) The same FOS spectropolarimetry and also the 
ground-based spectropolarimetry shows no evidence 
of any enhancement in polarization in the broad lines 
relative to the nuclear continuum. The broad lines 
are believed to originate only in a region much more 
compact than the scattering region that we have seen 
in our imaging polarimetry data. Any diluting radi- 
ation such as starlight or free-free emission from the 
scattering region itself (Tran 1995c), only dilutes the 
polarization of the scattered continuum, not the scat- 
tered broad lines. Therefore, if diluting radiation were 
present, the polarization would be higher in the broad 
lines than in the continuum, which has been observed 
in several Seyfert 2 galaxies, but not in NGC 1068 
(Antonucci et al. 1994; Tran 1995a,b). 

Those spectropolarimetric data, however, cover a 
large area, 4. "3 x l."4 aperture in Antonucci et al.'s 
data and a 2. "4 wide slit in Tran's data, which could 
mean that the assumption of scattered-light domina- 
tion might be true only in the brightest clouds. There- 
fore we restrict ourselves to use only the polarization 
degree of the bright knots. We will discuss the validity 
of this assumption again in 
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The following method for obtaining the three- 
dimensional gas distribution and its result have been 
discussed in Kishimoto (1997). We present the method 
and results here, with the errors in the observed polar- 
ization that have been estimated in much more detail, 
and with the position of the nucleus that has been ob- 
tained more accurately. 

6.2. Method 

The polarization of the scattered radiation from 
a cloud depends basically on its viewing angle, but 
also on the size of the cloud compared to the distance 
to the illuminating source. This effect is large if the 
cloud is very close to the source. In this section, we 
describe a simple method for interpreting the degree 
of polarization, taking this effect into account. 

Consider a spherical, uniform cloud of radius R, 
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which is optically thin to scattering. The cloud is il- 
luminated by a point source from a distance r (see 
Figure Ha). We consider only electron scattering, be- 
cause the primary scatterers in the innermost region 
of NGC 1068 are thought to be electrons (Antonucci 
& Miller 1985; Miller, Goodrich & Matthews 1991). 
The calculation of polarization for optically thin sys- 
tems of axisymmetric distribution has been imple- 
mented by Brown & McLean (1977). We follow their 
naming conventions for a geometrical factor a and 
shape factor 7. Our three quantities, polarization of 
the scattered radiation P, viewing angle ip, and the 
relative cloud radius 77 = R/r, are related as 



P = 



— sin ip 



2a (ry) -|- sin^ -0 ' 



where a{ri) is a function of 77 and defined as 

a=(l+7)/(l-37). 

We find the shape factor 7 as 

.,3 I _ ^2 



7 



2k{rf} 



where 



k{rj) — Tj ^ {1 — rj ) In ■ 



l + ?7 



(7) 



(8) 



(9) 



(10) 



Since the projected distance r^ = rsin-i/; is the ob- 
servational quantity whereas the actual distance r is 
not, we re- write -q as 



rj — r/sin'0, where 77' = R/rp. 



(11) 



These relations are rather complicated, but for small 77, 
7 becomes close to unity, hence a close to —1, and 
eq. (M) becomes the well-known equation for Thomson 
scattering as 



P = 



1 — COS^ ij} 

1 -1- cos^ -0 



(12) 



This means that the shape of the cloud or the distri- 
bution of the gas inside the cloud does not affect the 
polarization of the scattered radiation if the cloud size 
is small enough compared with the distance to the il- 
luminating source, and the polarization depends only 
on the viewing angle, as a natural consequence. 

The three-dimensional positions of the clouds can 
now be obtained, given the position of the central 
source in the image (xc, j/c)- If we denote the position 



of the center of each cloud in the image as (x, y), the 
distance z for each cloud from the sky plane in which 
the nucleus resides (see Fig. 0b) is written as 



where 



j/tani/j. 



Tp = \/{x-XcY + {y~VcY 



(13) 



(14) 



One of the ambiguities in this method is that, due 
to the forward-backward symmetry of the optically 
thin electron scattering, there are two possible view- 
ing angles derived from a single value of polarization, 
i.e. the forward and backward scattering case. These 
two angles correspond to the two positions along the 
line of sight, shown as Figure ^. Also, if a cloud 
is not sufficiently resolved, several clumps may exist 
along one line of sight. This could be the limit of 
our method, but in this case, the two positions de- 
rived from the polarization set a constraint on the 
positions of the clouds. If there is a scattering clump 
between these two positions, there must be another 
one outside of these two positions, to dilute the highly 
polarized scattered light from the inner clump. In this 
sense, these two positions should be considered as the 
"critical" positions. Having these ideas in mind, we 
simply assume that there is only one cloud along one 
line of sight. 

6.3. 3D mapping 

Given the position of the nucleus (xc, 27c) , the view- 
ing angle can be calculated numerically from the 
two observational quantities, polarization P and the 
relative size of cloud 77', using equations (|7|)^(|ll|). 
Figure |lO| shows the map of the viewing angle on the 
plane of (77', P) for a forward scattering case. There is 
an upper limit in P for certain 77' when ip is 90°. This 
maximum polarization becomes smaller for larger rj' . 
This is regarded as a "geometrical dilution" of po- 
larization, where the position angles of polarization 
of the scattered light from different places inside the 
cloud become significantly different so that the polar- 
ization tends to cancel out. As a consequence, there 
is a "region of uncertainty" at around 77' = 0.9 '--^ 1.0 
and P = 25 '--^ 35%, where two viewing angles can 
be solutions for a single point of (77', P). For the ob- 
servational points in this region, with certain values 
for the errors in 77' and P, we can obtain the viewing 
angle only with a large uncertainty. 

In Figure [l^, we have plotted the observational 
points of the clumps assuming that the cloud radius 
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Fig. 10. — The map of viewing angle ip on the plane 
of {rj',P). The observed positions of each cloud are also 
shown as squares with errors. 

is our aperture radius 0."15/2, with the total error in 
P estimated in § S, and the error in r]' that originates 
from the error in the position of the nucleus. We have 
taken the error of the nuclear position symmetrically 
as (-0."03±0."05, -0."12±0."10) for simplicity. Ex- 
cept the central few clouds, we see the viewing angles 
are well approximated by equation (12) which is the 
relation for r]' = 0. This means that the cloud shape 
or gas distribution inside the cloud does not affect 
the result as described above, although we have used 
equation (0) to derive the viewing angles for all the 
clumps. The polarization for each cloud is shown in 
Figure |ll|. The results are summarized in Table ^, 
with equatorial coordinates as in Table 0. 

The errors in ip and z are calculated only from the 
error in P in Table |3| and we did not include the error 
in the nucleus position, since it has a systematic effect 
on all clouds. The polarization degrees have been 
debiased following Simmons & Stewart (1985). That 
is, the final S/N in polarization, P/ap, where up is 
the total uncertainties in P as estimated in § 0, is large 
enough for all the clumps so that we are allowed to 
use the equation given in Wardle & Kronberg (1974), 
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Fig. 11. — Observed polarization degree of each cloud 
is shown on the center of each cloud in white numbers. 
The ID numbers for each cloud are also shown as black 
numbers in circle, displaced to the east and west for lower 
and upper positions, respectively. 



where Pobs is calculated simply as y^Q^ + U^/I. 

From these data, we have constructed the three- 
dimensional view of the nuclear region in Figure n2l 
Each cloud is illustrated as a sphere with 0."15 di- 
ameter. The position of the nucleus is shown as a 
smaller sphere with its error. The solid line at each 
cloud indicates the error of its position from the un- 
certainties in P. The positions in this figure would be 
only for one of the two possible cases. We could select 
another side of the critical positions for each cloud, 
or even both. However, as we also discuss in the next 
section, there is observational evidence that suggest 
an absorption excess in the southern region compared 
to the northern region, from a color difference (Mac- 
chetto et al. 1994) and HI absorption (Gallimore et 
al. 1994). Therefore we have simply selected the case 
shown in Figure ll^. 

7. Discussion 



(15) 
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Table 3 
Polarization and 3D position of cloud 



ID No. 


x(") 


y {"T 


R.A. (J2000) 


Dec. (J2000) 


P(%) 


epA(°) 


V>(°) 


zn 


1 


0.19 


-1.50 


02 42 40.698 


-00 00 49.31 


16.4 ±5.2 


96.6 ± 


8.9 


32.1 ±4.9 


2.23 ±0.42 


2 


0.23 


-1.34 


40.696 


49.15 


23.4 ±4.3 


104.9 ± 


6.6 


38.1 ±3.4 


1.60 ±0.19 


3 


0.54 


-1.26 


40.675 


49.07 


15.6 ±4.8 


122.3 ± 


9.0 


31.3 ±4.7 


2.09 ±0.38 


4 


0.36 


-1.20 


40.687 


49.01 


17.0 ±4.5 


103.5 ± 


8.5 


32.6 ±4.1 


1.80 ±0.29 


5 


0.19 


-1.18 


40.698 


48.99 


31.5 ±4.9 


95.4 ± 


5.1 


43.8 ±3.2 


1.12 ±0.13 


6 


0.28 


-0.87 


40.693 


48.68 


20.3 ±4.5 


106.4 ± 


7.4 


35.6 ±3.7 


1.14 ±0.16 


7 


-0.19 


-0.73 


40.724 


48.54 


43.7 ±5.9 


71.5 ± 


3.5 


51.4 ±3.4 


0.50 ±0.06 


8 


-0.19 


-0.49 


40.723 


48.30 


42.6 ±4.9 


62.8 ± 


3.8 


50.9 ±2.9 


0.32 ±0.03 


9 


0.00 


0.02'' 


40.711 


47.79 


19.8 ±4.3 


63.2 ± 


6.3 


36.0 ±3.9 


0.19 ±0.03 


10 


-0.07 


0.10 


40.716 


47.71 


17.3 ±4.8 


91.2 ± 


7.3 


33.2 ±4.5 


0.34 ±0.06 


11 


0.04 


0.24 


40.708 


47.57 


19.3 ±5.1 


86.4 ± 


6.4 


34.7 ±4.4 


0.53 ±0.09 


12 


-0.42 


0.30 


40.739 


47.51 


15.3 ±4.6 


138.4 ± 


7.5 


31.1 ±4.5 


0.96 ±0.17 


13 


-0.56 


0.31 


40.748 


47.50 


25.1 ±4.9 


139.2 ± 


4.9 


39.3 ±3.7 


0.84 ±0.11 


14 


0.11 


0.36 


40.704 


47.45 


19.3 ±5.3 


74.7 ± 


6.3 


34.7 ±4.5 


0.72 ±0.12 


15 


-0.12 


0.39 


40.719 


47.42 


13.6 ±4.2 


99.5 ± 


9.3 


29.3 ±4.3 


0.92 ±0.16 


16 


-0.29 


0.40 


40.731 


47.41 


11.6 ±4.4 


125.6 ±10.0 


27.2 ±5.0 


1.14 ±0.25 


17 


-0.58 


0.49 


40.749 


47.32 


19.9 ±4.8 


135.2 ± 


6.1 


35.2 ±4.1 


1.16 ±0.18 


18 


0.02 


0.52 


40.709 


47.29 


9.9 ±4.6 


80.3 ±11.0 


25.2 ±5.6 


1.37 ±0.35 


19 


-0.35 


0.55 


40.734 


47.26 


18.1 ±4.5 


121.7 ± 


7.2 


33.6 ±4.0 


1.11 ±0.17 


20 


-0.29 


0.65 


40.730 


47.16 


14.9 ±4.2 


111.4 ± 


9.0 


30.6 ±4.2 


1.37 ±0.23 


21 


-0.23 


0.77 


40.727 


47.04 


14.5 ±4.3 


111.9 ± 


9.5 


30.3 ±4.3 


1.57 ±0.27 


22 


-0.88 


0.81 


40.770 


47.00 


11.7 ±5.0 


127.0 ±10.9 


27.3 ±5.6 


2.44 ± 0.59 


23 


-0.39 


0.85 


40.737 


46.96 


20.0 ±4.1 


105.7 ± 


7.1 


35.3 ±3.5 


1.46 ±0.19 


24 


-0.58 


1.01 


40.750 


46.80 


42.4 ±3.8 


112.1 ± 


4.1 


50.5 ±2.2 


1.03 ±0.08 


25 


-0.13 


1.05 


40.720 


46.76 


19.3 ±4.9 


96.0 ± 


7.6 


34.6 ±4.2 


1.70 ±0.27 


26 


-0.73 


1.24 


40.759 


46.57 


16.9 ±4.3 


115.4 ± 


8.2 


32.6 ±4.0 


2.39 ±0.36 


27 


-0.80 


1.38 


40.765 


46.43 


23.2 ±4.2 


115.2 ± 


6.2 


37.9 ±3.2 


2.17 ±0.25 


28 


-0.59 


1.52 


40.750 


46.29 


27.2 ±4.1 


105.5 ± 


5.7 


40.9 ±2.9 


2.00 ±0.21 


29 


-0.68 


1.53 


40.756 


46.28 


26.0 ±4.1 


109.4 ± 


6.0 


40.0 ±3.0 


2.11 ±0.23 


30 


-0.19 


1.55 


40.723 


46.26 


19.9 ±5.3 


85.1 ± 


6.9 


35.2 ±4.5 


2.38 ± 0.40 



'^Measured from the UV peak, which is taken as the peak position in the 'sharp UV image (with fine focus). 

''The position of the maximum intensity in the / image (total intensity, out of focus) is slightly different from that 
in the 'sharp UV image'. 
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Fig. 12. — Three-dimensional view of the nuclear region. 
Each cloud is drawn as a sphere of diameter 0."15, and 
the nucleus is shown as a smaller sphere in the center, 
with the errors in its position. The scale unit is one arc- 
second. The projected views are also illustrated on each 
coordinate plane with each cloud drawn as an open circle. 
The direction to us is indicated by the large white arrow. 
The error in the position of the cloud center along the line 
of sight is indicated by the solid line for each cloud. 
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Fig. 13. — The relative UV flux for each cloud in Fig. O 
is plotted as its ID number, against the projected distance 
from the nucleus, in log scale. The horizontal lines repre- 
sent the error in the projected distance, where the solid 
(dotted) ones are for the case with the nucleus shifted to 
the north (south). The slope of —2 is indicated in the 
upper-right. 



7.1. PA map 

We have shown that the UV polarization image 
is completely consistent with a point-source scatter- 
ing within the accuracy of the FOC polarimetry, over 
wide regions from ~100 pc scale down to ~10 pc scale 
around the nucleus. The scattered radiation is very 
extended to the east and west, as well as to the north 
and south, as shown in Figure ||. We should note, 
however, that there could be much smaller deviations 
from a point-source scattering, which we cannot dis- 
cuss with the given limit of the accuracy of these 
imaging polarimetry data. We just certainly do not 
see any significant deviations from a centrosymmetric 
pattern beyond the estimated errors of the data, when 
we exclude the regions contaminated by the effect of 
the bad focus. 

Young et al. (1996) have estimated the size of the 
'torus' that is thought to be obscuring the nucleus. 
They have calculated its size to be greater than 200 pc 
(~ 3") in diameter, based on the absorption feature 
at ~ 1" south from the nucleus in their near-infrared 
polarized flux image. In Figure |l^, we have plotted 



the UV intensities of each cloud against the projected 
distance of each cloud from the nucleus, with the error 
of the projected distance originating from the error in 
the position of the nucleus. The statistical errors of 
the UV intensities were quite small, so are not shown 
in the figure. We clearly see that the fluxes of the 
southern clouds (Nos.l~8) are systematically smaller 
than those of the northern clouds, which suggests an 
absorption excess in the south. However, if we as- 
sume that the radial distribution of the UV flux is 
the same in the south as in the north, the observed 
UV flux ratio suggests an absorption of only Ay <^ 1 
at ~ 0."5 — 1" south from the nucleus. The obscuring 
'torus' could be extended to this scale, but it should 
be much smaller if we define it as the absorbing ma- 
terial which has enough optical thickness to hide the 
scattered radiation completely. The size of this ma- 
terial must be less than 0."3 projected onto the sky 
which corresponds to ~ 20 pc, based on the distribu- 
tion of the polarization seen in Figure ^. 

Muxlow et al. (1996) and Gallimore et al.(1996c) 
argue that the most probable location of the nucleus 
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7.2. 3D distribution 

The largest source of ambiguity in tlie derived 3D 
distribution is the assumption of scattered-hght domi- 
nation in the UV rang e. B oth of the two observational 
results described in § 6.1 are for large apertures and 



-0.5 0.0 

Offset West (arcsec) 



Fig. 14. — The registration of the 5-GHz MERLIN radio 
map (Muxlow et al. 1996) onto the HST shaxp UV image, 
using the astrometric result of Capetti et al. (1997). Both 
are in log scale. The registration uncertainty is shown as 
a dash-dot plus at the SI source. Our new location of the 
nucleus is drawn as a solid plus, whose size represents our 
estimated error. 

in the VLBI radio image is one of the most southern 
components, called the SI source (Gallimore et al. 
1996a). Strong water maser sources have been found 
at this SI source (Gallimore et al. 1996b; Green- 
hill et al. 1996; GreenhiU & Gwinn 1997), and these 
sources have been interpreted to be associated with 
the obscuring torus around the nucleus. In Table ^, 
we show the offset of the SI source from the UV peak, 
taken from the absolute astrometric result of Capetti 
et al. (1997). In Figure |l|, the 5-GHz MERLIN radio 
map by Muxlow et al.(1996) is registered on the sharp 
UV image (F253M filter, fine focus), using this astro- 
metric result. The uncertainty of the registration is 
indicated at the SI source. Our new location of the 
nucleus is shown with the error stated in Table ^. The 
new location coincides with the SI source within the 
errors. Therefore our result supports that the hidden 
nucleus resides in the SI source. 



our assumption may not be valid especially in the 
faint regions. 

We can roughly check this scattered-light domina- 
tion by looking at the radial dependence of the UV 
flux in Figure n3. If we assume that the gas density, 
volume, and viewing angle are the same for all the 
clumps, the plotted points should distribute with a 
slope of —2. The overall distribution mostly follows 
this slope. However, some of the clumps are above the 
line. This could be compensated by the differences 
in density and volume, but it also would suggest the 
existence of diluting radiation (non-scattered, unpo- 
larized light) in these clouds. 

The determination of the viewing angle of the UV 
brightest cloud (UV peak) could be beyond the limit 
of our method, since the cloud is too close to the nu- 
cleus. In this case, inhomogeneity of the cloud within 
the aperture will have a significant effect on the po- 
larization, while these imaging polarimetry data do 
not have enough spatial resolution for further inves- 
tigation. This uncertainty, however, does not affect 
the overall 3D distribution in Figure 11^. 

A correlation between the optical morphology in 
the HST images and the radio jet structure (Muxlow 
et al. 1996; Gallimore et al. 1996a, c; see Fig.|l4[) has 
been suggested by Capetti et al. (1997). Our 3D 
distribution obtained might indicate a certain linear 
structure, which could be related to the jet structure. 
However, much more detailed consideration of the na- 
ture of the UV radiation is needed. Detailed analysis 
using other HST images, including the discussion of 
the fraction of diluting radiation, will be presented in 
a subsequent paper (Kishimoto 1998). 

8. Conclusions 

We have analyzed the HST/FOC imaging polarime- 
try data of NGC 1068, and shown that the data 
are consistent with a simple point-source scattering 
model, within the FOG polarimetric accuracy. We 
have re-determined the location of the nucleus, by 
eliminating the regions which are suspected to have 
significant contaminations from the neighboring pix- 
els. The error circle suggests that the nucleus is lo- 
cated between clouds A and B, and the most probable 
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location of the nucleus has been found to be very close 
to cloud B, as close as ~ 0."08, which is much closer 
than has ever been claimed before. 

Based on this result, we have derived the three- 
dimensional view of the nuclear gas distribution, as- 
suming that the UV radiation is dominated by scat- 
tered radiation in the bright knots. The inferred 
three-dimensional distribution of the clouds might 
suggest the existence of a linear structure which could 
be related to the radio jet, though the assumption of 
scattered-light domination should be examined using 
other high-resolution images. 

The author thanks Robert Jedrzejewski for pro- 
viding various information on the FOC polarimetry, 
and Robert Antonucci for his helpful suggestions to 
improve the manuscript. The author would like to 
thank Ryuko Hirata for his kind advice and gener- 
ous support. The author also appreciates discussions 
with David Axon and Alessandro Capetti, and would 
like to thank Jack Gallimore for kindly providing the 
electronic data for the MERLIN radio map. 



A. Error Estimation 

In this appendix we summarize the method for the 
error estimation discussed in § |[ 

Consider the incident radiation of the Stokes pa- 
rameters {I,Q,U,V). The flux through three polar- 
izers, fi{i — 1, 2, 3), are written as 

f^ - ^^^iI - h COS 261, -Q^h sin 26*, • U), (Al) 

where ti,ki, and 9i are each polarizer's transmit- 
tance, polarization efficiency, and axis direction, re- 
spectively. We define the polarizers' axis directions to 
be measured counter-clockwise from the -\-x axis di- 
rection of the image, while we set the reference plane 
of the Stokes parameter along the y axis so that the 
position angle of the polarization is measured counter- 
clockwise from the y direction. 

The Stokes parameters (/, Q, U) are derived from 
the observed fluxes fi through each polarizer by the 
inverse relation of the above equation. For conve- 
nience, we define 



/; 



t, 



■h, 



(A2) 



which corresponds to the transmittance-corrected in- 
cident flux for unpolarized light. If we write (/, Q, U) 



also as {Ii, 12,13), these are calculated as 
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where 
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k2h sin(-26'2 + 263) 



k2 



Z2 sin 2^2 



Z2 sm z»2 -t- ^3 sin W3 

12 cos 2^2 ^ ^3 COS 26*3 

fcgfci sin(-26'3 4- 26*1) 
— fca sin 203 -I- k\ sin 2Q\ 



^3 sm AKJ3 -h K\ sm z(7i 
^3 cos 2^3 — fci cos 2Q\ 

- h2 sin(-20i + 2^2) 



k-yh 



n\H2 snu^— zai -t- ZP2^ 
—k\ sin 29 1 -t- fc2 sin 2^2 
k\ cos 29\ — ^2 cos 2^2 
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and 



A = 



fc2fc3sin(-26l2 + 26'3) 
-|-fc3fcisin(-26'3 + 26li) 
-|-fcifc2sin(-26'i + 26l2) 



(A5) 



Note that for the nominal case of 9\ = 180°, 6'2 = 
60°, 6*3 = 120°, and ki = k2 = ks = 1, this equation 
simply becomes 



(A6) 



For the FOC polarizers at the transmission range of 
F253M filter, we take the polarization efficiencies to 
be (Nota et al. 1996, p38). 
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fci = 0.986, fe = 0.976, fca = 0.973. 



(A7) 



The polarization degree P and position angle 0pA are 
calculated from I,Q,U as 



P = 



I 



1 u 

= — arctan -- . 

2 Q 



(A8) 



Now we estimate the total error in P with the fol- 
lowing four major sources; 



<^p = (c^p ) + \<^p ) +\<^p ) +['^p ) ■ (A9) 

The statistical error in the photon counts ap^^ is a 
random error, while the other three errors, each de- 
scribed below, would be systematic errors. We treat 
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these three, however, in the same way as for random 
errors, since we do not have enough caHbration results 
to correct for them exactly. 

The error from the image registration uncertain- 
ties, cr|J"^*, is estimated by actually shifting the images 
with the calibration uncertainties of 0.3 pixel (Hodge 
1995) and calculating the resulting change of P. The 



other errors arc calculated using equations (A2) 
(A8) in the following way. 



The statistical error ap'^* is written as 



^axis^ In particular, ct™® simply becomes 



1 






(A15) 



when we set the uncertainty in all the three polarizers' 
axes to be ag and all three polarization efhciencies to 
be the same. 
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where c/; is calculated assuming poisson noise in the 
counts. The error from the uncertainties in the direc- 
tions of polarizers' axes, a^^^, is calculated as 



^axis \ 2 



dP 



3 

E V do,. 



2 



(All) 



where we take ag^ to be 3° (Nota et. al 1996, p36; 
see also Robinson & Thomson 1995). Finally, afi'" 
is the error from the uncertainties in the "correction 
factors" , discussed in § 0, to be multiplied to each of 
fi. To take this into account, we re- write // to define 
correction factors S,i as 



Ji ~ vi ■ , Ji 



(A12) 



and derive dp ^^ from 



/_corr\2 



3 

E 

4=1 



dP 



(A13) 



where we take a^. at most to be 0.05 for the cases (A) 
and (C), 0.04 for the case (B), estimating from three 
kinds of uncertainties discussed in § y. A ll p art ial 
deriv ative s can be calculated from equations (A2)^(A8) 
and ( |A12|) , setting 6*1 == 180°, 02 = 60°, 6I3 = 120° and 
i^ = l. 

The error in ^pa can also be estimated in just the 
same manner as ap above. The well-known approxi- 
mate relation in high S/N case. 



cep 



1 (Tp 

2~P 



(in radian). 



(AM) 



holds only for a^p^'^ and a'^p" (note the similarity be- 
tween dP/df, and dP/d^i), 



but not for crfj"" 
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